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The transmission of a probe field experiencing electromagnetically induced transparency and opti- 
cal switching in an atomic medium enclosed in an optical cavity is investigated. Using a semiclassical 
input-output theory for the interaction between an ensemble of four-level atoms and three optical 
cavity fields coupled to the same spatial cavity mode, we derive the steady-state transmission spec- 
tra of the probe field and discuss the dynamics of the intracavity field buildup. The analytical and 
numerical results are in good agreement with recent experiments with ion Coulomb crystals [M. 
Albert et al, Nature Photon. 5, 633 (2011)]. 
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I. INTRODUCTION 

Electromagnetically Induced Transparency (EIT) is a 
quantum interference phenomenon occurring when two 
electromagnetic fields resonantly excite two different 
transitions sharing a common state [H-Q- intense 
control field addressing one of the transitions can sub- 
stantially modify the linear dispersion and absorption of 
an atomic medium for a weak probe field resonant with 
the second transition. Since its first observation by Boiler 
et al. [3] , EIT has been successfully exploited for instance 
to control the propagation of light pulses through an oth- 
erwise opaque medium for light storage and retrieval ^ 
[T2j and quantum memories [l3l - l22| . Besides providing 
a means for controlling the linear susceptibility of an 
atomic medium EIT can also be exploited for generat- 
ing strong optical nonlinearities 0, H, [H, . For in- 
stance, in the four-level atomic configuration such as the 
one depicted in Fig. [T]d, the nonlinear susceptibility of 
the medium can be strongly enhanced at the same time 
as the linear susceptibility is suppressed. The large cross- 
Kerr effect between the probe field and a third switching 
field can then be used e.g. for high-efficiency photon 
counting psj . all-optical switchi ng f26l[27l and nonlinear 
optics at low-light levels 0, H, uM. l28l |29|. nonclassical 
state generation [3^ or the realization of strongly inter- 
acting photon gases [3l|, • 

When such a nonlinear EIT medium is positioned in 
an optical cavity one first-of-all benefits from the en- 
hanced interaction of the ensemble with well-defined 
spatio-temporal field modes. This is of great value for 
enhancing the effective optical depth and for realizing 
high-efficiency quantum memories [33l - l36| or Fock state 
quantum filters The EIT- induced reduction of the 
cavity linewidth |38| can also be used to increase the sen- 
sitivity of atomic magnetometers t2Ll39|, enhance cavity 
optomechanical cooling processes p0l44a . achieve lasing 
without inversion [43 [jOptical switching |44l. |45| or quan- 
tum state swapping 46j . 

On the other hand, the cavity EIT-interaction can be 



furthermore exploited to enhance the cross-Kerr non- 
linearity to investigate e.g. photon-blockade mecha- 
nisms |47H49| . photon-photon interactions [50], highly- 
entangled state generation (sij or novel quantum phase 
transitions for light [s^ . 

Cavity EIT has been observed with atomic beams 53 , 
in cold and room temperature atomic ensembles (s^ - 56| , 
and even with single or few atoms in high-finesse optical 
cavities [13, [13 • Cavity EIT as well as EIT-based optical 
switching have also recently been observed with cold ion 
Coulomb crystals [soj . 

We theoretically investigate here the interaction be- 
tween an ensemble of atoms with the four-level structure 
depicted in Fig.[T)D and three optical fields coupled to the 
same spatial cavity mode. Using a semiclassical input- 
output theory we derive the atomic linear and nonlinear 
susceptibilities for the probe field and its steady state 
transmission spectra as it experiences EIT or an EIT- 
based cross-Kerr effect. We study in particular the effect 
of the transverse mode profiles of the fields on the normal 
mode spectrum of the atom-cavity system as in the ex- 
periments of Ref. [soj and compare with the more usual 
situation where only the probe field is coupled to the 
cavity 

The paper is organized as follows: in Sec. |TT] we in- 
troduce the system under consideration and derive the 
equations of motion. We obtain the probe field suscep- 
tibility and transmission spectrum, first when the probe 
and control fields interact with the atoms in an EIT situ- 
ation in Sec, mil then when all three cavity fields interact 
simultaneously with the atomic medium in Sec. IIVI Fi- 
nally, all-optical switching at low-light levels using ion 
Coulomb crystals as in the experiments of Ref. [59| is 
discussed in Sec. IVl 



II. SYSTEM CONSIDERED AND EQUATIONS 
OF MOTION 

We consider an ensemble of Ntot four-level atoms with 
the level structure depicted in Fig. [I] where levels |1) 
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FIG. 1. (a) Atomic four-level structure considered, (b) Sys- 
tem considered: an ensemble of four-level atoms, positioned 
in a linear optical cavity, interact with three optical fields. In 
the all-cavity geometry all three fields are injected into the 
cavity and coupled to the same cavity mode. In an alterna- 
tive standard geometry, only the probe field ap is injected into 
the cavity, while the control and switching fields, Oc and as, 
axe free-propagating with waists much larger than that of the 
ensemble (see text for details). 

and 1 2) are long-lived ground or metastable states, and 
levels |3) and |4) are excited states. The atoms are en- 
closed in a linear optical cavity where they interact with 
three optical fields: a probe field on the |1) — > |3) tran- 
sition, a control field on the |2) — > |3) transition and a 
switching field on the |2) — > |4) transition. The cavity 
is assumed to have asymmetric mirror transmissions, Th 
and Tl, with Th ^ Tl, and, unless stated otherwise, 
all fields are injected into the cavity through the high- 
transmission mirror. The annihilation operators of the 
field modes in the Heisenberg picture and in the rotating 
frame are denoted by dp, dc and ds, respectively, with 

the convention [aa,aj,] = I, {a = p,c,s). a\ii, = 
is the atomic operator associated with the j— th ion po- 
sitioned at Tj {^,v = 1 — 4). We assume that all three 
fields are resonant or close to a resonance with the same 
spatial cavity mode, which we take to be the fundamen- 
tal Gaussian TEMqo mode for simplicity. The interaction 
Hamiltonian in the rotating-wave approximation and in 
the rotating frame is given by 

^<^f = ~ ^ XI 5p*p(ri) cos [kpZj ^ iPp{Yj)] ap4i^ (1) 
j 

- ?i^5c*c(rj) cos [kcZj + ipc{rj)] dca^£ 

3 

- cos [ksZj + ipsirj)] d^a^j] + h.c, 
j 

where the g's are the maximal single-atom coupling 
strengths for the transitions considered, the ^''s and (p'a 
are the fields' transverse mode functions and longitudi- 
nal mode phases, respectively [60j . Since our main goal 
is to discuss the effects related to the spatial transverse 
structure of the fields we will consider for simplicity the 
situation where the length, 2L, of the ensemble is much 
smaller than the Rayleigh range of the cavity, as e.g. in 
the experiments of Ref. [s^. We thus neglect the lon- 
gitudinal variations of the phases and the waists of the 
light fields over the length of the ensemble. We there- 



fore set the longitudinal mode phases to and assume 
a Gaussian transverse structure of the fields given by 
"i'ai'r) = exp(— r^/w^) (a = p, c, s), where Wa is the cav- 
ity waist considered. We will in the following examine 
two situations: 

(i) The novel case in which all three fields have the 
same transverse mode profile and the ensemble has a 
large radial extension as compared to the waists, like e.g. 
in the experiments of [5^ . We shall refer to this situation 
as the all-cavity case. 

(ii) The more usual situation in which the control and 
switching fields have a large transverse intensity profile as 
compared to the extension of the ensemble. This would 
typically be the case if these fields were interacting with 
the atoms not through the cavity (53 - l58| . or for an en- 
semble radially confined to a region with dimension much 
smaller than the waists, as could be obtained with e.g. 
with a string of atoms or two-component ion Coulomb 
crystals ^611]. We shall refer to such a situation, in which 
the transverse mode profiles of the fields can be ignored, 
as the standard case. 

We will in addition assume that, because of their mo- 
tion, the atoms are "warm" enough such that they probe 
any field variation along the longitudinal standing-wave 
structures of the fields during the characteristic time 
scales of the dynamics of the fields due their interactions 
with the atoms. As discussed e.g. in [62| - [65j and in 
Appendix |X1 one can under these conditions assume av- 
eraged longitudinal couplings (ja — ga/V^ [a = p,c,s). 
Keeping only the transverse spatial dependence of the 
cavity modes, the Hamiltonian becomes 

j 

+ gs'^sirj)d,sa'fj+h.c. (2) 

For comparison the "cold" atom situation where the 
atoms are well-localized with respect to the field stand- 
ing wave structures during the interaction is treated in 
Appendix [X] 

The atom-field dynamics of the observable mean values 
can be standardly derived via 6 — {l/ih){[d,H]), where 
o = (6) is the mean value of observable 6 and the total 
Hamiltonian H = Ha + Hf -I- Haf is the sum of the 
interaction Hamiltonian ^ and of the atomic and field 
Hamiltonians 

Ha = -hJ^Apag^ + {Ap - A,)ag^ + {Ap - A, + A,)aii\ 
j 

Hf = -hApdldp - hA^dldc - hA^^dlds, (3) 

where Ap = ojp — ui^i , Ac — — (^32 and Ag = cj^ — W42 
are the one-photon detunings, and A^ ~ cua ~ i^cav are 
the cavity detunings between the fields with frequency 
u!a and the cavity resonance frequency cocav (a = p, c, s). 
Denoting by 731, 732 and 742 the spontaneous decay 
rates and introducing a phenomenological decay rate 
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7o for the ground-state coherence operators (tJj (70 ^ 731,732,742), one obtains the fohowing set of coupled dif- 
ferential equations 



^12 = -(70 - j(Ap - Ac))CTr2 - Wp^p{rj)ap4i> + igc^cirj)al<7\y + ig,^ s{rj)a*a\'^ (4) 

(^13 = -(7 - i\)<^i3 + «5p*p('^j>p(^n - + igc'i'c{rj)aca[i^ (5) 

= "(7s - «(Ap + As))ctJ^4^ - i.9p*p(rj)apCT^^4^ + i.9,*,(rj>^cri^2^ (6) 

0-23 = -(7 - i^c)(^23 + *5c*c(?'j)ac(CT^2^ - a^^^) + igp'i'p{rj)apai\^ - igsilJs{rj)asa^^^ (7) 

0-24 = -(7s - i^s)(^2l - igc'^c{rj)ac(j''^} + igs"^ s{rj)as{a''^^ - cr^^l) (8) 

0-34 = -(7 + ls- 270 - »(As - ^c))<ji'I - igc-^c{rj)ala'^^} + ig,*,(rj>,4^2' - Wp'^p{rj)a*p<^il (9) 

= 731^3' - *5p*p(^j)ap4l' + *5p*p('^j>;f^?3' (10) 

0-22^ = 732^33' + 742cri-'2 - igc^c{rj)ac(j''^^ + ig^-^ rirj)alcT^^^ - ig^-^ s{rj)as(J^^^ + ig^-^ s{rj)a*(T^^} (11) 

<^33^ = -(731 + 732)cr33^ + igp^ p{r j)apai^} - igp^p{rj)a*(T[^^ + igc-^c{rj)accr^£ - igc*c(rj)a*cr^^3^ (12) 

(jii = -7i2crii + igs'i'sirj)as(7^^ - ig s{r j)ala^il (13) 

ap = -(« - zApap -riY^-gp^p(r,)a^l -I- ^^aj," (14) 

a, = -{K- iA^)a, + i ^ gc*c(r,)ag) + ^^a^ (15) 

d. = -(«: - *A^)a, + z ^ -gs^s{r,)a'il + ^^ar (16) 



where 7 = (731 -I- 732)72 -I- 70, 7s = 742/2 -|- 70 and r is 
the cavity round-trip time. The input fields are denoted 
by (a — p,c,s). The total cavity field decay rate 
(assumed equal for all fields for simplicity) is denoted by 
K = kh + kt + ka, where kh,l = Th,l/'^t are the de- 
cay rates corresponding to the mirrors' transmission and 
= ^/2t is a decay rate corresponding to round-trip 
absorption losses A. While these absorption losses are 
not essential to the understanding of the physical mech- 
anisms studied here, we include them for completeness 
as they ver y o ften affect experiments with high-finesse 
cavities I57H59I. 



The previous set of equations can be solved numerically 
for any initial internal atomic state, input field pulses, 
ensemble geometry and atomic distribution. We focus in 
the following on the situation in which all the atoms are 
in state |1) initially and the probe field is much weaker 
than the control and switching fields, so that one can 
perform a first-order expansion in the probe field to get 
analytical expressions for various quantities, such as the 
probe susceptibility, its cavity transmission and reflec- 
tion, the EIT buildup time, etc. 



III. CAVITY ELECTROMAGNETICALLY 
INDUCED TRANSPARENCY 

A. EIT regime 

We first investigate the EIT situation where the atoms 
interact with both the probe and control fields, but no 
switching field is injected into the cavity. The input 
probe and control fields are abruptly switched on at time 
t = and have thereafter constant intensities. We place 
ourselves in the weak probe regime, when gp|ap| ^ ffclocl 
and the intracavity photon number is much smaller than 
the number of interacting atoms. All the atoms are then 
essentially in |1), and the only non-zero atomic compo- 
nents at first order are the probe optical dipole a^^ and 
the ground-state coherence 0] ■ We also assume that 
the control field is tuned to resonance with the |2) — > |3) 
transition (Ac — 0) and the cavity is resonant with the 
|2) — !■ |3) transition, i.e. Ap = A^ = (5 = A. As 
the control field probes no atom, its intracavity ampli- 
tude reaches its steady state value in a time . Since 
we are interested in getting simple analytical expressions 
for the steady state of the system and its dynamics over 
the typically slower EIT buildup tiniescales, we can con- 
sider that the control field intracavity Rabi frequency is 
constant and equal to its steady state value fic = gcO-c- 



4 



Equivalently, the control field can be turned on slightly 
before the probe pulse is applied. The relevant equations 
of motion governing the evolution of the intracavity probe 
field are then 

dp = -{k ~ iA)ap + igp'^'^p{rj)a[^^ + ^J2^^a}^ , 



(17) 



'l3 



2 ' 



U) 



(18) 
(19) 



B. Steady state susceptibility 



These equations can be readily solved in steady state 
to obtain the mean value of the intracavity probe field, 



iA 



iXEIT 



where the EIT susceptibility is given by 



Xeit 



7- 



iA- 



70 -i A 



(20) 



(21) 



For large ensembles of atoms with non-correlated posi- 
tions, we can reformulate Eq. ((2T|) in terms of the local 
atomic density p{r). Assuming the same transverse pro- 
files for the control and probe fields ^'p(?') — 'i'c{r) — 
exp(— r^/w^), Eq. (|2ip can then be recasted in 



^9 



XEIT = — ^ / dr p(r) 



-2r-'/w^ 



7- 



iA- 



.(22) 



where flc = gc\o-c\- For an ensemble with a uniform 
atomic density, such as ion Coulomb crystals in a lin- 
ear radiofrequency trap as used in [H, [66l - l68| , and with 
a large radial extension as compared to the cavity waist 
(see e.g. [59, 66-68]), the integral can be calculated an- 
alytically and yields a susceptibility 



Xeit 



iglN ln(l + e) 



where 



e 



iA 



e 



(7--^A)(7o - iA) 



(23) 



(24) 



is an effective saturation parameter for the two-photon 

2 

transition and = P^^L is the effective number of 
atoms defined in [g^, [68l |. This result can be compared 
to the standard situation in which the control field waist 
is much larger than that of the probe field [ij, H, 111] 



St _ ^9pN 
Xeit 



1 



7 - iA 1 -h e 



(25) 
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FIG. 2. (Color online) Real (blue) and imaginary (red) 
parts of the probe first-order susceptibility as a function of 
the probe detuning A (solid line: all-cavity EIT [Eq. ((23])], 
dashed line: standard EIT [Eq. 1)25^ ]'). The parameters, 
{gpVN, 7, 70, He) = 2nx (16, 11.2, 6 x 10^", 6)MHz, are sim- 
ilar to those used in the experiments of Ref. [591 ]. 



The real and imaginary part of these susceptibilities 
are plotted in Fig. [5] for typical parameters used in the ex- 
periments with ion Coulomb crystals of [5^. They show 
the typical transparency window in the absorption pro- 
file and the rapid change in dispersion around two-photon 
resonance. In comparison with the standard one, the all- 
cavity susceptibility clearly shows non-Lorentzian line- 
shapes, as expected from its different dependence with 
respect to Q [Eqs. and ((25t ]. 

From Eq. (PH]) and the input-output relations 

re f 



V2 



v 



(26) 



the steady state cavity transmission and reflection for the 
probe field are given by 

.tr 2 ^ ,— — 2 



p 



R: 



a: 



p 



kh + i^L + HA — iA — ix 
kh - kl - KA + iA + ix 



kh + hl + ka ~ iA — ix 



(27) 
(28) 



Using Eqs. (123125^ one can then compute the normal 
mode spectrum of the probe field transmission. In the 
collective strong coupling regime, when gp\/N > k,7, 
one expects three normal modes in the transmission spec- 
trum: two modes at probe detunings ±JgpN + Vl1/2, 



corresponding to the two-level ensemble modes it^p V N 
shifted by the presence of the control field, and one mode 
at zero-detuning for the probe (two-photon resonance 
here), corresponding to the cavity EIT resonance [55, 69']. 

This is illustrated in Fig. |31 where the probe trans- 
mission spectra, normalized to the bare cavity resonant 
value 

2 



Tn = 



iA 



IX 



(29) 
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FIG. 3. Upper panel: Probe field normalized transmission 
spectra for an empty cavity (dotted line), a cavity containing 
a uniform density ensemble interacting with the probe field 
only (dashed grey line) and a cavity containing a uniform 
density ensemble interacting with a probe and a control field 
in an all-cavity (solid black line) and standard EIT situation 
(gray solid/dashed lines). Parameters: ((jp^iV, 7, 70, Sic, «:) = 
27rx(16, 11.2, 6x10"", 6, 2.2) MHz. The solid and dashed grey 
curves show the standard EIT situation with control field Rabi 
frequencies ilc ~ {2tt)6 MHz and fic = iic,eB = (27r)6/2.2 
MHz, respectively. Lower panel: Spectra enlarged around 
A = 0. As discussed in the text, the latter value of the ef- 
fective Rabi frequency was chosen to illustrate the lineshape 
difference in a situation when the all-cavity and standard 
EIT resonance curves have comparable halfwidths at half- 
maximum. 



is represented for the cases of (i) an empty cavity (x = 0), 
(ii) Bi cavity containing a uniform density ensemble inter- 
acting with the probe field only (x = WpN/ (7 — «A)) and 
a cavity containing a uniform density ensemble interact- 
ing with a probe and a control field in an all- cavity (Hi) 
and standard (iv) EIT situation, for which the probe sus- 
ceptibility is given by (f23| and (|25|) , respectively. One ob- 
serves indeed that the value of the cavity transmission in 
presence of EIT is restored to close to the bare cavity res- 
onant value in a narrow frequency window around reso- 
nance. The width of the central EIT feature can be calcu- 
lated by expanding the transmission around two-photon 



resonance. In the standard case and in the regime con- 
sidered previously (770 <^fll< gpN), one finds that the 
probe transmission is Lorentzian-shaped around A = 
(see e.g. [i,|3|) 



To: 



9lN 



1 



a 



\keit - iAI 



(30) 



The interaction thus emulates a cavity with an effective 
halfwidth 



K-EIT 



7o 



(31) 



which is smaller than the bare cavity halfwidth k when 
gp^/N > fie- The analysis is a bit more complicated in 
the all-cavity case, due to the complex dependence of the 
susceptibility with the saturation parameter G and the 
non-Lorcntzian profile (as can be seen e.g. from Fig.[3]D). 

The different dependence of the susceptibility with 
the effective saturation parameter Q makes it in general 
impossible to define an effective control field Rabi fre- 
quency in the all- cavity situation which would give the 
same susceptibility or transmission as in the standard 
case. However, if one is interested in comparing situa- 
tions in which the EIT resonance features have similar 
widths, one can perform a similar expansion of the nor- 
malized transmission given by Eqs. (|23l) and p9| around 
A = and define an EIT resonance width also in the 
all-cavity situation. In the regime 51^ 3> 770, the all- 
cavity EIT resonance width matches the standard EIT 
one with an effective control field Rabi frequency ilc,cff — 
na/[{y'2TT^ + 4(ln2C)2-7r)/(7rV2-H2(ln2C)2)]i/2^ ^here 
C — glf/2Kj. For the parameters of Fig. O the scal- 
ing factor for the effective Rabi frequency is ^ 2.2 for 
instance. On the other hand, if one was interested in 
comparing the minimum absorption level on two-photon 
resonance , one could define an effective control Rabi fre- 
quency as n^gg ~ r2c/\/lii(ric/77o); in the EIT regime 
where 770 and for a given value of fie. The reso- 

nant absorption in the all-cavity situation with a con- 
trol field Rabi frequency fie could then be effectively 
compared to that of a standard situation in which the 
maximal Rabi frequency has been scaled by a factor 
^ v/hr(fi2/77o). 

Fig.|4]shows the corresponding reflectivity spectra for a 
cavity having non negligible roundtrip absorption losses, 
as observed in the experiments of (59.] . In general, since 
the reflected fleld results from the interference between 
the input fleld and the intracavity fleld, the reflectiv- 
ity levels have a slightly more complex dependence on 
the atomic absorption and the cavity losses. The reflec- 
tivity spectrum exhibits nonetheless the same qualita- 
tive features as the transmission, with two normal modes 



at frequencies g^N -\- and a third one at zero- 
two-photon detuning corresponding to the reduction of 
atomic absorption due to the EIT effect. Effective con- 
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and we have assumed the initial conditions ap{0) 
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FIG. 4. Reflection spectra for the same configurations and 
parameters as in Fig. [S] for a cavity with non negligible 
roundtrip absorption losses [{kh , i^l, ha) = (27r) x (1.5, 0, 0.7) 
MHz]. 



trol field Rabi frequencies can also be defined in a similar 
fashion as previously, as shown e.g. in [sol [70j. 



C. Dynamics 

In this section we focus on the dynamics towards 
reaching the steady state during a resonant EIT inter- 
action (A = 0). Assuming again a constant control field 
Rabi frequency and performing a Laplace transform of 
Eqs. P?)l . (HH]) and yields the following equations 



(k + s)ap[s] = i 5p*p(^j)c^i3 [s] 



2kh 



where the Laplace transform of / is defined by 



(J)r 



f[s] = £[f{t)] = / dt. 

Jo 



'fit) 



(32) 



(7^-^^(0) = o-\y{0) = 0. These equations allow for ex- 
tracting the Laplace transform of the intracavity probe 
field amplitude 



ap[s] 



E 



(33) 



j 7+s+|f^cP/(7o+s) 



and calculate its time evolution by performing the inverse 
Laplace transform. It is however instructing to look at 
the dynamics in the adiabatic limit in which the effective 
cavity linewidth emulated by the EIT medium is smaller 
than the bare cavity linewidth and the dipole decay rate, 
i.e. Kbit < ^,7. In this limit it can be shown that the 
intracavity field and the optical coherence adiabatically 
both follow the ground state coherence, which evolves 
at a rate keit- In the standard case, from Eq. p3p . 
one finds that the intracavity field amplitude increases 
exponentially with a time constant 1/keit, consistently 
with the steady state spectrum analysis of the previous 
section. In the all-cavity case the inverse Laplace trans- 
form has to be calculated numerically. It yields a non- 
exponential increase in the intracavity field intensity oc- 
curring on a timescale approximately given by 1/keit, 
with Qc scaled as in the previous section. 




10 15 20 

Time (^s) 

FIG. 5. Probe field normalized transmission as a function of 
time (solid black line: all-cavity case, solid grey line: standard 
case, dotted grey line: empty cavity) for a resonant probe 
field (A = 0). Parameters as in Fig.O The control field Rabi 
frequency is scaled by a factor 2.2 in the standard situation. 

Figure [S] shows the time evolution of the normalized 
probe transmission in the two situations discussed, for 
the same parameters as in the previous section and for 
an input probe pulse abruptly switched on at i = and a 
constant control field. For comparison the (much faster) 
bare cavity response is also shown. 



IV. OPTICAL SWITCHING 

We now turn to the all-optical switching situation, in 
which the transition |2) i — > |4) is addressed by the 
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switching field a^, while the control and probe fields are 
in an EIT situation. When the switching field is detuned 
from atomic resonance (|As| ^ 75) and weak enough 
such that the absorption to level |4) is negligible, its main 
effect is to light-shift level |2), thereby changing the bare 
EIT resonance condition for the control and probe fields. 
When the light-shift becomes comparable or greater than 
the width of the cavity EIT window, the transmission of 
the probe field is inhibited, as the cavity is switched off 
resonance by the presence of the switching field. 



A. Probe susceptibility 




-0.4 -0.2 0.2 0.4 

Probe detuning A/In (MHz) 

FIG. 6. (Color online) Optical switching suscepti- 
bility: real (blue) and imaginary (red) parts around 
two-photon resonance. The solid line shows the all- 
cavity switching for ((jpx/iV, 7, k, 70, f2c, 7s, As, fis) — 2tt x 
(16,11.2,2.2,0.0006,4,11,4300,40) MHz. The dashed line 
shows the standard switching situation for the same parame- 
ters, but with the control and switching field Rabi frequencies 
scaled by a factor 2.2 for comparison. 

We assume again that almost all the atoms stay in 
|1) and that the control and switching field intracavity 
Rabi frequencies have reached their steady state values, 
ric = (jcCic and rig — (jsO-s, when the probe is injected. 
Performing a first-order treatment in the probe field, the 
equations of motion for the non-zero coherences are given 
by 



■0) 

'^13 



(i) 



(35) 

(3) 



(34) 



"14 



(36) 



Solving Eqs. ([M)) . (1551) and (pS)) in steady state readily 
yields a mean intracavity probe field amplitude of the 
form (PU)) . with a susceptibility 



Xsw = 2j*fp*p(''j) 



70 



iA- 



7s — iAs — 



(37) 



with = (7s|as|- For a uniform density medium with 
an extension larger than the cavity waist and for fields 
with identical transverse profiles, one gets an analytical 
expression for the susceptibility 



Xsw 



where 



7 — iA 



e. = 



ein(i + e + e,s) 



(e + e. 



e + e. 



{-fs - iAs - iA){"fo - iA) 



(38) 



(39) 



is defined analogously to the effective EIT saturation pa- 
rameter for the probe. This susceptibility can again be 
compared to that of the standard case where the control 
and switching fields have waists much larger than that of 
the probe 



1 



7-iAi + e/(i + e,)' 



(40) 



B. Probe field transmission spectrum 

The real and imaginary part of both susceptibilities 
are shown in Fig. |5] for typical parameter values taken 
from (59| . As expected, in the standard case, the ef- 
fect of the switching field is to shift the position of the 
EIT resonance by an amount that corresponds to the 
AC Stark shift of level |2). The effect is more com- 
plex in the all-cavity case, as the shift for each atom de- 
pends on its radial position, thus leading to an asymmet- 
ric frequency behavior of the absorption and dispersion 
around two-photon resonance. This effects are manifest 
on the probe field transmission spectra, which are shown 
in Fig. [7] for different switching field intensities. While, in 
the standard configuration, the probe transmission pro- 
file is shifted away from the bare two-photon resonance 
without too much distortion as the switching field inten- 
sity increases, it is substantially distorted in the all-cavity 
case due to the different AC Stark shifts experienced by 
atoms at different radial positions. The accuracy of these 
analytical expressions for the susceptibility and transmis- 
sion have been checked by numerically solving Eqs. 
([TS]) . These findings are also in good agreement with the 
experimental observations of Ref . [5^ . 



V. LOW-LIGHT LEVEL OPTICAL SWITCHING 
WITH ION COULOMB CRYSTALS IN CAVITIES 

We now turn to the prospects of achieving low-light 
optical switching of a single-photon probe field using ion 
Coulomb crystals in an optical cavity and base our prac- 
tical discussion on the parameters of Refs. [H, H^. We 
assume an asymmetric linear cavity geometry similar to 
that described in [s^, with length ^ 1 cm and finesse ~ 
4000 and n ~ hh = (27r)1.5 MHz. We consider an inter- 
action with "'"Ca+ ions on the id '^I?3/2, wj = -1-3/2 — 
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photon number is then calculated in steady state, yield- 
ing the cavity transmission. Using a control field Rabi 
frequency Q.c = (27r)2 MHz allows for increasing the res- 
onant probe transmission to ^ 90%. The variation of 
the probe transmission for different switching field de- 
tunings As and intracavity photon numbers Ug is shown 
in Fig. [8] under these conditions. For all these simula- 
tions we checked that the absorption of photons to level 
1 4) was negligible and that the depletion of atoms from 
level |1) remained at most at the percent level. Optical 
switching is observed to take place with increasing pho- 
ton numbers as the detuning is increased, as expected 
from the previous discussion and analysis. We define 



! 0.8 



i 0.6 



T3 0.4 
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FIG. 7. (Color online) Probe field normalized transmission 
spectrum for different switching field Rabi frequencies Q.s : (a) 
All-cavity switching, fi^ = 27rx (0, 15, 25, 40) MHz from left to 
right. Other parameters as in Fig. [6l (b) Standard switching 
for the same parameters but with Rabi frequencies Q,c and Q,s 
scaled by a factor 2.2. 



Ap 2Pi/2,mj = +1/2 (probe), 3d 3^^3/2, mj = -1/2 ^ 
4p ^A/2,wj = +1/2 (control) and id ^ZJg/a,™,/ = 
— 1/2 — >• Ap ^P3/2,wj = -1-1/2 (switching) transitions, 
for which the respective maximal ion coupling strength 
are {gp,gc,gs) = (27r) x (0.53,0.22,0.18) MHz respec- 
tively. We assume a standard situation and numerically 
calculate the steady state normalized probe transmission 
by solving Eqs. P|)- ([T5|) for a probe input field intensity 
such that the mean intracavity photon number is one in 
steady state in an empty resonant cavity. 

Taking an effective collective coupling strength g\/N = 
(27r)16 MHz renders the crystal/cavity system com- 
pletely opaque for the probe field in the absence of 
control field (To ~ 1%). We assume that the cavity 
field decay rate is the same for all fields and that the 
atomic decay rates as defined in (|4l)- (fT6l) are (7,7s,7o) = 
(27r) X (11.2, 11, 6 X lO"'^) MHz. For the simulations the 
control and switching fields are injected 0.5 ^s before the 
probe field, with rise times much shorter than the inverse 
of the cavity field decay rate, to allow them for reaching 
their steady state values. The mean intracavity probe 




1 10 100 1000 10" 

Mean number of intracavity switching photons 



10' 



FIG. 8. (Color online) Probe field normalized transmission as 
a function of intracavity switching photon number Us for dif- 
ferent switching field detunings A3 (blue: A3 = (27r)0 MHz, 
red: A3 = (27r)110 MHz, green: A3 = (27r)4300 MHz). Pa- 
rameters: (gp^iV, 7, 73,7o,«:, r^c. A) = 27r X (16, 11.2, 11,6 x 
10"'', 1.5, 2,0) MHz. 



the minimal switching photon number n* as the mini- 
mal number of intracavity switching photons needed to 
bring the normalized transmission from 90 % to 10 %. 
In the standard situation, one can easily show from (j25p 
and (|29| that having a 90 % transmission in FIT im- 
poses that keit ^ 4O70. To get substantial switching 
we require that the light-shift induced by the switching 
field rig /2 As is a few times the width of the FIT trans- 
parency window Kbit- A numerical estimation shows 
that f2s/2As ~ 5keit, which gives n* ~ AOOjoAg/g'^. In 
agreement with Fig. [5] we find that ~ 17000 photons are 
needed for the large detuning of 4.3 GHz used in [5^ and 
~ 400 for a detuning of ~ IO7S. The previous estimate is 
actually still valid for a resonant switching field replacing 
As by 7s, which would give a minimal photon number 
of ^ 40 for the parameters of Fig. [S] This illustrative 
numerical example is based on the experimental param- 
eters of [59| . but we note that lower switching numbers 
could in principle be reached e.g. using smaller cavities 
or stronger switching transitions 
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VI. CONCLUSION 

Using a scmiclassical theory for the interaction of 
four-level atoms with three optical cavity fields, the 
effect of the transverse mode profiles on the susceptibility 
and transmission spectrum of a probe field experiencing 
EIT or EIT-based optical switching has been discussed. 
Contrarily to the standard situation where the control 
and switching field Rabi frequencies are the same for 
all atoms, non-Lorentzian EIT resonance lineshapes and 
asymmetrical switching lineshapes are predicted when all 
three fields are coupled to the same cavity mode. Closed 
analytical forms for the susceptibility and transmission 
spectrum of the probe field have been found to explain 
these lineshapes, in good agreement with numerical 
simulations and with experiments using ion Coulomb 
crystals in cavities. Last, the prospect for achieving 
low-light optical switching with ion Coulomb crystals in 
moderate finesse optical cavities was discussed. 
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FIG. 9. (Color online) Probe field transmission spectrum 
around two-photon resonance for the three situations: stan- 
dard EIT (dotted line) and all-cavity EIT with localized 
(dashed line) and delocalized (solid line) atoms. Parameters 
as in Fig. [3] with the control field Rabi frequency scaled by a 
factor 2.2 in the standard situation. 



one can recast (jAip into 



"31,+ + "31,- 



(A2) 



h.c. 



Appendix A: "Localized" or "delocalized" atom 
situations 



In this section we investigate the effect of the atomic 
motion on the cavity EIT feature. Free-space EIT 
with standing- wave field geometries has been investigated 
both theoretically [gl, [gI, [H, [71| - [73| and experimen- 
tally [H, [131 , and the influence of the atomic motion 
on the storage and retrieval of pulses in such geometries 
has been discussed in e.g. |63 - l65| . 

To discuss the influence of the atoms' longitudinal ve- 
locities on the cavity fleld spectrum in an EIT situation 
in a relatively simple fashion, one can start from the 
Hamiltonian ([T]) without switching field and neglecting 
the variation of the control and probe field Gouy phases 
and fc- vector difference (fcj, ~ fcc = /c) 



Haf = - h gp^pjrj) cos(fczj-)apO-^{^ 
j 

— fey~^ge^c(i'j) cos{kzj)acd-^^2 +h.c.. 



(Al) 



Introducing ± 



(t\^^ exp(±ifczj) 



and (7 



23, ± 



(T23 exp(±ikzj) the optical dipole operators correspond- 
ing to the two running wave fields propagating with ±k, 



Making the same assumptions as in Sec. IIIIl i.e. almost 
all atoms in state |I), constant and strong control field 
Rabi frequency, one gets equations of motion in an EIT 
situation which are similar to Eqs. (|I7II8II9p . but now 
include the longitudinal dependence of the coupling with 
the fields: 



do) _ 

'^13, ± — 



■A3) 
'12 



(A3) 

(7 " ^A)ag]^ + Z5p*p(r,)ap(f + e^'^''^')/2 

(A4) 



(70 - ^A)a^ + ^^^:^cir,)iaY,'.+ + o\'l_)/2. 

(A5) 



(j) , A]) 



If the typical timescales for the longitudinal atomic mo- 
tion (trapping frequencies, thermal motion,...) are faster 
than the EIT dynamics timescale, but still slower as com- 
pared to the atomic dipole dynamics (i.e. if the typical 
longitudinal velocity v is such that keit ^ l^^^l ^ 7), 
then the terms in exp(±2ifcz) can be averaged out in (IA4p 
and one retrieves the " delocalized" situation discussed in 
this paper or in the experiments of [59| . Physically, this 
can be explained by the fact that the moving atoms will 
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only be in two-photon resonance with the co-propagating 
parts of the standing waves. As the atoms hence see on 
average fields with a longitudinal intensity which is half 
the maximum of the standing wave value, their dipole is 
reduced which means that the coupling strengths can be 
effectively rescaled by l/\/2, yielding the effective Hamil- 
tonian 

If we now assume that the atoms are sufficiently cold 
for their longitudinal positions to be fixed with respect 
to the longitudinal standing-wave structure of the cavity 
fields during the EIT interaction, one keeps the longi- 
tudinal dependance in the coupling terms to solve the 
previous equations of motion, which are equivalent to 



-{k — iA)ap + igp ^p(^j) cos{kzj)a 



13 



Solving Eqs. (jA6IA7IA8|) in steady state yields a suscep- 
tibility 



_ ng'Pe(r-3)^cos^(fcz,-) 



(A9) 



For a large, uniform density ensemble with random lon- 
gitudinal ion positions along the cavity axis, one gets 



Xeit 



ig;N 21n[(l-H Vl-H29)/2] 
7 - iA e 



(AID) 



'^13 



'12 



(A6) 



(7 — iA)a[^:] + igp^p{rj) cos{kzj)ap 

(A7) 



-(70 - «A)crp2^ + ini-^drj) cos{kzj)a[^^ . 



(A8) 



where 8 is given by Eq. (12^ . The probe field trans- 
mission spectrum around two-photon resonance is shown 
in Fig. [S] for the same parameters as in Fig. ^jp in the 
three situations considered: standard EIT and local- 
ized/delocalized all-cavity EIT. The localized situation 
is seen to give rise to a slightly broader EIT resonance, 
since the atoms see on average a slightly higher effec- 
tive control field Rabi frequency, as one would intuitively 
expects. 
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